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An array feed combining system for the recovery of signal-to-noise ratio (SNR) 
loss due to antenna reflector deformation has been implemented and is currently be- 
ing evaluated on the Jet Propulsion Laboratory 34-m DSS-13 antenna. The current 
signal-combining system operates under the assumption that the white Gaussian 
noise processes in the received signals from different array elements are mutually 
uncorrelated. However, experimental data at DSS 13 indicate that these noise pro- 
cesses are indeed mutually correlated. The objective of this work is to develop a 
signal-combining system optimized to account for the mutual correlations between 
these noise processes. The set of optimum combining weight coefficients that maxi- 
mizes the combined signal SNR in the correlated noises environment is determined. 
These optimum weights depend on unknown signal and noise covariance parameters. 
A maximum-likelihood approach is developed to estimate these unknown parame- 
ters to obtain estimates of the optimum weight coefficients based on residual carrier 
signal samples. The actual combined signal SNR using the estimated weight coef- 
ficients is derived and shown to converge to the maximum achievable SNR as the 
number of signal samples increases. These results are also verified by simulation. 
A numerical example shows a significant improvement in SNR performance can be 
obtained, especially when the amount of correlation increases. 


I. Introduction 

An array feed-combining system has been proposed for the recovery of signal-to-noise ratio (SNR) 
losses caused by large antenna reflector deformations at Ka-band (32-GHz) frequencies in the Deep 
Space Network [1]. In this system, a focal plane feed array is used to collect the defocused signal fields 
that result from these deformations. All the signal power captured by the feed array is then recovered 
using real-time signal-processing and signal-combining techniques. The optimum combiner weights that 
maximize the combined signal SNR were derived in [1] under the assumption that the white Gaussian 
noise processes in the received signals from different array elements are mutually uncorrelated. These 
optimum weights depend on unknown signal and noise parameters that need to be estimated. The work in 
[1] proposed to estimate the optimum weights from the observed residual carrier received-signal samples 
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using a maximum-likelihood (ML) approach. The actual combined-signal SNR in an uncorrelated noises 
environment when the estimated weights are used in place of the optimum weight coefficients was also 
derived in [1]. 

A seven-element array feed combiner system is currently being evaluated at the JPL DSS-13 34-m 
antenna. Although the work in [1] assumed mutually uncorrelated noise processes, experimental data [2] 
indicate that the noise processes in the received signals from different feed elements are indeed correlated, 
with correlation coefficients of the order of 0.01 under clear sky conditions. Since the noise in each of 
the array feed element signals consists of receiver white noise plus noise due to background radiation, it 
has been conjectured that this small correlation is caused by near-field atmospheric background noise. 
Recent data gathered at DSS 13 in more adverse weather conditions, however, indicate both increases 
and decreases in the observed amount of correlation. This experimental work is still in progress. Larger 
correlations may also result from undesired radiation source emissions gathered by the antenna side lobes. 
Moreover, our recent work [5] has shown that the array feed combining system derived in [1], which is 
suboptimal in the correlated noises environment, actually can have a better performance in the presence 
of correlated noises. Therefore, it is important to develop optimum signal-combining techniques that 
account for the mutual correlation between the noises in the signals from different array elements. That 
is the objective of this work. 

As a first step towards this objective, in Section II we provide a derivation of the set of optimum 
combining weight coefficients that maximizes the combined-signal SNR in this correlated noises environ- 
ment. These optimum weights depend on unknown signal amplitude and phase parameters as well as 
noise variance and correlation parameters. An ML approach is then developed in Section III to estimate 
these unknown parameters and arrive at an ML estimate of the optimum weight coefficients based on 
residual carrier received-signal samples. The actual combined-signal SNR is derived in Section IV when 
the estimated weights are used in place of the optimum weights, with the details given in Appendices A 
and B. The SNR performance is shown to converge to the maximum achievable SNR as the number of 
signal samples used in the estimates increases. These results are also verified by simulation. Numerical 
examples are given in Section V with a particular choice of noise covariance matrix and signal parame- 
ters and show a significant improvement in SNR performance compared to the previous signal-combining 
system developed in [1], especially when the amount of correlation increases. 


II. Array Feed Signals and Optimal Combining Weights 

Consider a A-element array and the NASA Deep Space Network standard residual carrier modulation 
with binary phase shift keyed (PSK) modulated square- wave subcarrier [4]. The received signal from each 
array element is downconverted to baseband and sampled. Similar to the combining system proposed in 
[1], only the residual carrier portion of the received signal spectrum will be used to estimate the unknown 
parameters in the combiner weights. The full-spectrum modulated signals from the array elements, which 
contain both the modulated sidebands as well as the residual carrier spectrum, are subsequently combined. 
As in [1], assume that the higher bandwidth primitive baseband signal samples are lowpass filtered by 
averaging successive blocks of M B samples to yield a full-spectrum signal stream B for each array element. 
Additive white Gaussian noise is assumed to be present in the primitive baseband signal sequences from 
each of the array elements. The white Gaussian noises in the primitive baseband samples from different 
array elements are assumed to be mutually correlated. Specifically, the noise samples corresponding to 
different array elements are assumed to be mutually correlated at any given time instant, but uncorrelated 
at different time instants. Let 

Vk{iB) = V k [cos 6 + js(i B ) sin 6} +n k (i B ), i B = 1,2,- (1) 

denote the stream B signal samples from the fcth array element. The complex signal parameters Vfe, 1 < 
k < K, represent the unknown signal amplitude and phase parameters induced by the antenna reflector 
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deformation. Moreover, 5 is the modulation index, s(i B ) = ±1 is the transmitted data, and {n k (i B )} 
is the zero-mean white Gaussian noise corruption in the stream B signal samples from the fcth array 
element. The primitive baseband signal samples are also more narrowly lowpass filtered by averaging 
successive blocks of Ma samples to yield a residual carrier signal stream A for each array element. 
Clearly, Ma > Mb and r? = Ma/Mb is the ratio of the bandwidth of stream B to stream A. Let 

u k (i A ) =V k cos 6 + m k (i A ), ^a — 1,2,--- (2) 

denote the stream A signal samples from the A;th array element. Here is the zero-mean white 

Gaussian noise corruption in the stream A signal samples from the fcth array element. Let A*, A 1 ' , 
and A t denote the complex conjugate, the transpose, and the complex conjugate transpose of the 
matrix A, respectively. In order to specify the correlations between the white noise sequences cor- 
responding to different array elements, consider the noise vectors n(i B ) = (ni(i B ) • • • nx(i B )) T and 
UBJAa) — (^i(*a) • • • Then {n(i B )} and {to(m)} are each sequences of independent identi- 

cally distributed (i.i.d.) zero-mean complex Gaussian random vectors of dimension K. The respective 
covariance matrices R B = {r B kj} = E[n(is)n(is) t ] and R A = {rAkj} = E[m(i / i)m(* J 4 ) t ] of n(i B ) and 
m(iA ) specify the mutual correlations between the white noises in the signal streams from different array 
elements. For example, r Bk j is the correlation between the noise variables n k (i B ) and rij(i B ). Because 
of the different averaging rates in streams A and B on the primitive baseband signals, it follows that 
R b = t\R a - Finally, these different averaging rates also imply that m(M) is independent of n(i B ) pro- 
vided that ia < i-B and the samples averaged to yield m(iA ) occurred prior to the samples averaged to 
yield n(i B ). 

Application of the complex combining- weight coefficients i W k , 1 < k < K yields the combiner output 
sequence z c (i B ) — s c (i B ) + n c (i B ), where 


K 

Sc (i B ) = J2 w ^ eMiB)s 

k = 1 


(3) 


and 


K 

n c (i B ) = 'Y^W k n k {i B ) 
k = 1 


(4) 


are the signal and noise components, respectively. Define W = ( W\ ■ ■ ■ Wk) f to be the vector of 
combining-weight coefficients. The objective is to determine the optimum weight vector W that maxi- 
mizes the SNR of the combiner output defined by 


, W \ = l E Mfg )]| 2 lsc(fi ?)| 2 
Var[z c (i B )] Var[n c (i B )] 


(5) 


The optimum weight vector and the maximum achievable SNR have been derived previously in [3]. For 
the sake of completeness, we provide a derivation below that is slightly different from that in [3] . Define 
V — (Vj • • • Vk) T to be the vector of complex signal parameters. Then, from Eqs. (3) and (4), we have 
|s c (*b )| 2 = \W T V\ 2 and Var[n c (i B )] = W t R f W*. substitution of which in Eq. (5) yields 


7(00 = 


\W T V\ 2 

W t R b W* 


( 6 ) 
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Since R B is a positive definite Hermitian matrix, there is a unitary matrix Q such that 

Rb = $D 2 Q (7) 

where D is a real-valued K x K diagonal matrix with the fcth diagonal term given by y/r B kk- Using 
Eq. (7), we have 


W t R b W* = (w T Q^Dj ( DQW *) 

= (DQW*y (DQW*) = \\DQW *\\ 2 


Moreover, since Q is the inverse of the matrix , 


(m*) 


T' 


-1 



= D~ l Q 


Hence, by using Eq. (9), we can write 


W T V = W T (DQ*f \(DQ*y 




= (DQW*y (D~ l QV) 


( 8 ) 


(9) 


( 10 ) 


So substituting Eqs. (8) and (10) in Eq. (6) and applying the Schwartz inequality gets 


7(10 = 


I (DQW*y (D-'QV) I 2 
\\DQW *\\ 2 


WDQW-Tm^QVW 2 

\\DQW *\\ 2 


= \\Dr l QV \\ 2 =1MAX 


( 11 ) 


Moreover, equality holds in Eq. (11) if and only if for some complex-valued constant a, DQW* = 
a (D~ 1 QVy . So, by using Eq. (7), the set of optimum weight vectors W OPT that achieves 7 (W OPT ) = 
Imax is given by 


Wqpt = a ( &DT 2 Q ) V - a (R* B )~ l V* = ^ CSk) -1 1* (12) 

where a is an arbitrary complex-valued constant. Moreover, it follows from Eqs. (7) and (11) that the 
optimum SNR 7 max can be written as 
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1MAX = (D-'QVy (D-'QV) = V?R b 1 V 


(13) 


Note that, in the uncorrelated noises case, R B l — D 2 . So the set of optimum weight vectors, Eq. (12), 
in this case is given by W OPT — aD_~ 2 V*, and the optimum SNR is given by 

7 MAX = WD-'VW 2 = E {r^ = 7 (14) 

tZ i rBkk 

which is the sum of the array element output SNRs. These results for the uncorrelated noises case agree 
with previous results derived in [1]. 


III. Parameter Estimation 

The signal parameter vector V_ and the noise covariance matrix R A are not known and need to be 
estimated to obtain an estimate of one of the optimum weight vectors, W OPT . given by Eq. (12). Assuming 
that these unknown parameters are not random, we propose to use ML estimates based on the stream A 
residual carrier-signal vector samples {u(iji)}, where u(z^) = (u\(i A ) ■ ■ ■ UK(iA)) T ■ Instead of estimating 
V directly, consider estimating X = VcosS. Note from Eq. (2) that (u('m)} is an i.i.d. sequence 
of complex Gaussian random vectors with mean X and covariance matrix R A . It then follows from 
multivariate statistical analysis [6,7] that, based on observations {u(iA — 1), • • • ,u(i A — L)}, 

x ML {i A ) = \ £ a(0 ( 15 ) 

l=iA ~ L 

and 

^A,Mi(^) = ^ E [n(0 — ^Ml(^)] [n(0 _ =^Ml(^)] (16) 

1—lA — -k 

are the respective ML estimates of X and R A . By the invariant property of ML estimators, 

Y-ml^a) = ^-^Kml^a) (17) 

is then the ML estimate of V and 

W-opt^a) = — (E-aml^a)} Vml^a) (18) 

is the ML estimate of W OPT given by Eq. (12). 

It is well known [6,7] that X_ Mpj (i A ) is a complex Gaussian random vector with mean X_ and covariance 
matrix (1 /L)R A . Furthermore, X_ml(^a) and R mb (ia) are statistically independent, and LR ML {i J \) has 
the same distribution as the random matrix 
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( 19 ) 


i= 1 

where Z_ i is a sequence of i.i.d. zero-mean complex Gaussian random vectors with covariance matrix R A . 
This type of distribution is called a complex Wishart distribution with parameters R A and (L — 1), and A 
in Eq. (19) is said to have a CW(R A , L — 1) distribution [7]. It has been shown in [8] that, for L > K + 1, 

E W-‘I - « 


for a j K x K CW(R A ,L — 1) distributed random matrix A. Since (A * ) ~ 1 = (A -1 )*, it then follows from 
using the property of Eq. (20) for the complex Wishart matrix LR ML {i A ) along with Eqs. (12), (17), and 
(18) that, for L > K + 1, 


E 








olL 


V(L-K-I) 


R A ~ l V* 


L 

L-K-l 


W-OPT 


(21) 


Hence, the ML estimate W OPT (i , if of the optimum weight vector W OPT is actually a biased estimate. 
This is not a problem, since it is clear from Eq. (12) that the optimum weight vectors are not unique and 
any complex scaled version of an optimum weight vector is also optimum. So the constant a in Eq. (18) 
can be set arbitrarily. We shall set a = (L — K — l)/L for the purpose of normalization and also assume 
that L > K + 1. The ML estimate of the optimum weight vector that will be used here is, therefore, 


with mean 


Kml^a) 


L-K-l 

Lr) 

L-K-l 
Lr] cos 6 


(&A ,ml (U))~ Vul^a) 
(J^a,ml (M))" 





( 22 ) 


(23) 


that is also an optimum weight vector. 

Note that V_ ml (ia) and Ra,m l(m) are both consistent estimates, i.e., they both converge with prob- 
ability one to their respective expected values V and R A in the limit as L tends to infinity. So it follows 
that W_ ML (i A ) also converges with probability one to the optimum weight vector W Q in the limit as the 
number L of samples tends to infinity. The analysis in Appendix C shows that this convergence also holds 
in the mean-square sense. These properties indicate that we may expect the actual combiner output SNR 
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using the estimated weight vector W ML (i a) to " also converge to the maximum achievable SNR 7 max as 
L tends to infinity. That result will be shown in the next section, which derives an explicit expression for 
the actual combiner output SNR for finite L. 

As in [ 1 ], these weight coefficient estimates are used in a sliding- window structure to produce the 
following combiner output sequence: 


„ T 

z c {iB) = W ML (i A )y(i B ) 


(24) 


where i A is the largest integer less than i B , so that the residual carrier-signal vector samples 
{u (i A — l) , • • ■ , u (i A — L) } used for estimating W ML ( i A ) occur before the full-spectrum signal vector 
sample y(i B ) = {yi{i B ), • ■ • , y/c(*e)) • This ensures that the noise vector n(i B ) in y(i B ) is statistically 
independent of {u (i A — l) , ■ • • ,u(i A — L) }, and hence the statistical independence between W ML (i^) 
and n(i B ). 


IV. SNR Performance Analysis 

The combiner output, Eq. (24), can be written as 


z c (* b ) = s c (i B ) + n c (i B ) 


(25) 


where 


s c (i B ) = W T ML (i A )Ve^ s 


(26) 


and 


n c {i B ) = W_ ML (i A ) n(i B ) 


(27) 


are the signal and noise components, respectively. In the following analysis, V , R B , and s(i B ) are assumed 
to be nonrandom parameters. As noted above, W ML (i A ) and n(i B ) are statistically independent. Since 
the components of n(i B ) are all of zero mean, it follows from Eqs. (26) and (27) that n c (i B ) also has zero 
mean and moreover is uncorrelated with s c (i B ). Thus, it follows from Eq. (25) that the actual SNR 7 ml 
of the combiner output can be written as 


7 ML 


1 E[*c (is)] I 
Var[z(i B )] 


1 e M^b)] I 

Var[s c (j B )] + Var[n c (iB)] 


Now, it follows from Eqs. (26) and (23) that 


E[s c M] | 2 = | E [W T ML (i A )]V | 2 = | WlV | 2 = | V}R~ b x V I 2 = r MAX 


(28) 


(29) 


where 7 max is the maximum achievable SNR given by Eq. (13). The explicit expression, Eq. (A-9), 
for Var[s c (i B )] is derived in Appendix A. Moreover, the derivation in Appendix B yields the expression, 
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Eq. (B-5), for Var [n c (i B )}. So, by using Eqs. (29), (A-9), and (B-5) in Eq. (28), one arrives at the 
following expression for the actual combiner output SNR performance: 


7 ML 


7 MAX 

1MAX{1 + Cl) + C-i 


(30) 


where 


LK — K 2 — K + 1 (L- A- !)(£-!) f 1 \ 1 

1 (L - K)(L - K - 2) (L - K)(L -K- 2) V^cos 2 6 ) L-K- 2 1MAX 

K(L-l)(L-K-l) 1 

2 L(L - K){L - K - 2) r]cos 2 6 K 1 

Since both C\ and Ci converge to zero as L tends to infinity, it follows from Eq. (30) that the actual 
combiner output SNR 7 ml converges to the maximum possible achievable SNR, 7 max, as the number 
of signal samples, L, used in the estimates tends to infinity. 


V. Numerical Example 

We consider here a numerical example using a K = 7 element array feed. In this example, a modulation 
index, S = 80 deg, and a primitive sample period of To = 2.5 x 10~ 8 s are assumed. The full-spectrum 
modulation signal is assumed to be of bandwidth 2 x 10 6 Hz, which yields M B = 20. Moreover, the 
ratio of the full-spectrum bandwidth to the residual-carrier bandwidth is 77 = Ma/Mb = 200. A nominal 
P T /N 0 of 65 dB-Hz is considered with a corresponding 7 = ( Pt/N 0 )MbT 0 (recall that 7, which is given by 
Eq. (14), is the sum of the array element SNRs). The white Gaussian noise processes in the received signals 
from different array feed elements are assumed to be correlated. Moreover, the correlation magnitudes 
are assumed to be inversely proportional to the distances between feed centers (this assumption has not 
been verified for the array feed system at DSS 13). In the numerical examples below, the following noise 
covariance matrix is considered: 

H B = 


- 

1 

p e A 4 7 r + | £) 

t'max'' 

p e -i(^+4>) 

p e -^T+« 

t'max ^ 

p e -i(¥+« 

t'max'-' 


Pma*e 

Pm a 

e j( i f+4>) 

1 

P e -j( s r+4>) 

t'max 

£max, p~j ( Of- +(f>) 

y/3 

Pmax „j(2y-—4>) 
2 e 

£max_p~j(^f-+<l>) 

■n/3 

/We'A^+tfO 

Pmax 

e j(^+4>) 

p e A¥+« 

t'max'-' 

1 


Pmax 

vx 

Bass. 

2 e 

BmS^. pj{^T ~4 > ) 
y/S 

p max e j< - B ? +4,) 

vt e 

n 

t'lriax 

1 

P 

t'max'-' 

Pmax g 

v3 

pmax p~ 
2 e 

Ana,eA¥+« 

Bass e -j( 2 f--<l>) 
2 e 

pmax 

V3 6 

p g — 

t'max'-' 

1 

p e -A^+« 

t'max'-' 


Pmax 

e i(^ 2L +« 

fjnax. pj ( S f-+<t>) 

f>max 
2 e 

Pmax gj ( +(J>) 

y/3 


1 

A 'maxe j( -^~^ 

\p eA 8 ” 

[f'max v 

p e A¥+« 

Pmax p—j 

V3 6 

Pmax r l(2Z-+<t>) 

2 e 

Pmax pj ( +</>) 

VI e 

p g — — <£) 

f-'max'-' 

1 


(32) 


Figure 1 shows the K = 7 element feed array geometry and the relative distances between feed centers. 
The main feed is labeled feed 1 and is surrounded by the six others. In the correlation matrix, Eq. (32), 
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Fig. 1. Feed array geometry 
for K = 7 feeds. 


the noise power in the feed-element stream B received signal samples are all assumed to be equal and 
normalized to one ( r B kk — 1 f° r all k). The maximum possible noise correlation (and also correlation 
coefficient) magnitude is denoted by p max ■ The corresponding noise correlation magnitudes in this matrix 
are identical and reflect an inverse dependence on distances between feed centers. The noise correlation 
magnitudes between nearest neighbor feed pairs are equal to the maximum correlation magnitude p max 
(for example, between the main feed, feed 1, and any of the outer feeds). As can be seen from Fig. 1, 
the next nearest neighbor feed pairs (for example, feed 2 and feed 4) are of a distance equal to \/3 times 
the distance between nearest neighbor feed pairs, and hence have a noise correlation magnitude equal to 
Pmax/VS- Finally, the furthest feed pairs (for example, feed 2 and feed 5) are of a distance equal to twice 
the distance between closest feed pairs, and so have a noise correlation magnitude equal to p max /2. The 
parameter <j> in Eq. (32) specifies the noise correlation phases. The rationale for assigning these phases 
will be described below. Since the noise power in each feed element is equal to one, the sum of the feed 
element SNRs is equal to the total received power, and so 7 = Pt- For the complex signal parameters 
vector, F, we shall use 


y : 


a/ 


' 6 


■j&l L 

e J * 





6 


e 3 



(33) 


where 0 < (3 < 1 represents the fraction of total received signal power in the main antenna feed (feed 1) 
as defined in [9]. The remaining total received signal power is then evenly distributed among the other 
six feed elements (this assumption does not appear to be generally valid for severely distorted antenna 
reflectors). The signal phases in Eq. (33) were chosen arbitrarily. The parameter <p in Eq. (32) specifies 
the relation between the noise correlation phases and the signal parameter phases as follows: Let Vk 
denote the fcth component of the vector V given by Eq. (33). Then for k > j, <j> = [phase of rskj] 

— [(phase of Vk) — (phase of V))]. Moreover, since R B is a Hemitian matrix, —<j> = [phase of r B kj] 

— [(phase of Vk) — (phase of V})] when j > k. That is, (j) represents the difference between the phase of 
the noise correlations and the phase differences between the corresponding signal components in V. The 
significance of this will be discussed below. 

In the numerical examples below, the signal combining scheme given by Eq. (24) will be referred 
to as the correlated noises algorithm. The signal-combining scheme from [1] will, on the other hand, be 
referred to as the uncorrelated noises algorithm. A comparison of the SNR performance of these two signal 
combining schemes will be made below, assuming the noise covariance matrix given by Eq. (32) and the 
signal parameter vector given by Eq. (33). Our previous work in [5] determined the SNR performance of 
the uncorrelated noises algorithm in the environment where the noise processes in the feed element received 
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signals are indeed correlated. For a given signal vector V and set of noise correlation magnitudes, the 
results in [5] show that the best-case performance of the uncorrelated noises algorithm occurs when each 
of the noise correlations between feed pairs has a phase that is exactly 180 deg from the phase difference 
between the corresponding signal components in V. For the correlation matrix, Eq. (32), this “best case” 
scenario for the uncorrelated noises algorithm occurs when the parameter (p — 180 deg. Moreover, the 
worst-case performance was shown in [5] to occur when each of the noise correlations between feed pairs 
has a phase that is exactly equal to the phase difference between the corresponding signal components in 
V. This “worst-case” scenario for the uncorrelated noises algorithm occurs when the parameter cj> = 0 deg 
for the correlation matrix, Eq. (32). The numerical examples below will consider these two extreme cases 
only. The performance comparisons between these two algorithms can then be made under both the 
most-favorable and the least-favorable situations for the uncorrelated noises algorithm. 

The SNR performances will be compared in terms of the combining gain, which is the ratio of the 
actual SNR performance to the sum of the SNRs of the array feed element received signals ( 7 ml /7 
for the correlated noises algorithm). Note that 7 is also the maximum achievable SNR in the uncor- 
related noises environment. Hence, the combining gain represents the SNR gain relative to the best 
possible performance in the uncorrelated noises environment. We shall, therefore, refer to an SNR 
gain when the combining gain is positive (in dB) and to an SNR loss otherwise. The SNR perfor- 
mance analysis given in Section IV showed that the combining gain 7ml/7 of the correlated noises 
algorithm converges to the maximum possible combining gain, 7 max /l, in the limit as the number of 
samples, L , approaches infinity. Figure 2 plots this maximum possible combining gain as the phase 
parameter (p in Eq. (32) is varied between 0 and 180 deg for values of p max equal to 0.1, 0.15, and 
0.2, and when (3 = 0.7. These results show that the best- and worst-case scenarios for the correlated 
noises algorithm are the same as those for the uncorrelated noises algorithm. These results show that 



Fig. 2. Maximum possible combining gain versus with /? = 0.7 
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the maximum possible combining gain is positive for a majority of the <j> phase values. Moreover, the 
maximum possible combining gains increase with increasing p max . Simulations were also performed to 
validate the analytical result, Eq. (30), that yields the combining gain for the correlated noises algorithm. 
Table 1 compares the simulated combining gain to the analytical result for the best-case R B matrix (with 
<t> — 180 deg) at various values of p max and L. Table 2 displays the corresponding comparisons for the 
worst-case R B (with cf> — 0 deg). The simulated combining gain from these two tables can be seen to be 
within 3 percent of the analytical results. This appears to validate the analysis. 


Table 1. Correlated algorithm simulated 
combining gain for best-case R B . 


SNR Gain, dB 

L 

Pmax 

Simulation 

Analytical 

500 

0.050 

0.85896 

0.80340 

500 

0.100 

2.24048 

2.17891 

500 

0.150 

4.75655 

4.67706 

500 

0.175 

7.25759 

7.14633 

500 

0.200 

15.06133 

14.60054 

5000 

0.050 

0.88286 

0.87982 

5000 

0.100 

2.26456 

2.25971 

5000 

0.150 

4.78108 

4.77260 

5000 

0.175 

7.28267 

7.26961 

5000 

0.200 

15.08713 

15.03570 

10 6 

0.050 

0.86770 

0.88821 

10 6 

0.100 

2.24674 

2.26859 

10 6 

0.150 

4.76168 

4.78312 

10 6 

0.175 

7.26309 

7.28324 

10 6 

0.200 

15.06514 

15.08578 


In Fig. 3, the combining gains of both algorithms are plotted versus the number of samples, L , for 
(3 = 0.7 (70 percent of the power in the main feed), using the best-case R B and with p mox having values of 
0.05, 0.1, and 0.15. We notice that for a fixed pmax » the correlated noises algorithm (CNA in the figures) 
actually has a smaller combining gain than the uncorrelated noises algorithm (UNA in the figures) for 
small values of L below a threshold value. In the examples considered, this threshold value increases 
with decreasing p max and is always less than about L — 500. Note that although the ML estimator is 
asymptotically optimal as L — » oo, it is not necessarily the best estimator for small values of L. Hence, 
the correlated noises algorithm may not have the best possible performance at small values of L. The 
convergence of the combining gain to the optimum theoretically achievable value was proved in the last 
section. These examples show that convergence to within 0.1 dB of the limiting value occurs at about 
L = 1000 samples. Figure 4 shows the respective combining gains as a function of f3 , the fraction of total 
received signal power in the main feed, for L = 5000 samples and values of pmax equal to 0.05, 0.1, 0.15, 
and 0.2. We can see that the performance of the correlated noises algorithm is much less sensitive to (3 in 
this example than that of the uncorrelated noises algorithm. The performance superiority of the corre- 
lated noises algorithm also increases significantly with increasing (3. The respective combining gains versus 
Pmax with L = 5000 and (3 taking on values of 0.7, 0.8, and 0.9 are shown in Fig. 5. This figure shows that 
the performance improvement of the correlated noises algorithm over the uncorrelated noises algorithm 
increases significantly with increasing p max - The results in these three figures indicate that there can be 
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Table 2. Correlated algorithm simulated 
combining gain for worst-case fig. 


SNR Gain, dB 

L 

Pmax 

Simulation 

Analytical 

500 

0.050 

-0.63659 

-0.68802 

500 

0.100 

-1.05467 

-1.10542 

500 

0.150 

-1.33373 

-1.38419 

500 

0.175 

-1.43041 

-1.48083 

500 

0.200 

-1.50179 

-1.55220 

5000 

0.050 

-0.61306 

-0.61410 

5000 

0.100 

-1.03140 

-1.03182 

5000 

0.150 

-1.31079 

-1.31072 

5000 

0.175 

-1.40765 

-1.40738 

5000 

0.200 

-1.47920 

-1.47877 

10 6 

0.050 

-0.62206 

-0.60598 

10 6 

0.100 

-1.03724 

-1.02373 

10 6 

0.150 

-1.31348 

-1.30265 

10® 

0.175 

-1.40878 

-1.39931 

10 6 

0.200 

-1.47881 

-1.47070 
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Fig. 5. Combining gain versus p max with L = 5000 for best-case Jfe. 
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a substantial SNR gain in this best-case scenario, particularly for the correlated noises algorithm at larger 
values of p mox - We also note that the value of p ma x cannot exceed 0.2 in order to preserve the positive 
definiteness of R B . 

Figures 6 through 8 repeat Figs. 3 through 5, respectively, using the worst-case scenario R B matrix in- 
stead, while keeping the other parameters unchanged. These examples show that this worst-case scenario 
can result in an SNR loss, and the loss nominally increases with increasing p max . Figure 6 again shows 
that convergence of the combining gain for the correlated noises algorithm to within 0.1 dB of its limiting 
value occurs at about L = 1000 samples. The small sample performance of the correlated algorithm 
again lags that of the uncorrelated algorithm. Similarly to Fig. 4, Fig. 7 shows that the performance 
superiority of the correlated noises algorithm increases with increasing (3. However, as /3 approaches one, 
the SNR loss of the correlated noises algorithm turns around and starts to decrease with increasing p max . 
At large values of (3 close to one, the correlated noises algorithm can in fact have an SNR gain even 
when the uncorrelated noises algorithm still has an SNR loss. Finally, Fig. 8 shows that the SNR loss 
of the correlated noises algorithm in this unfavorable situation deteriorates much slower than that of the 
uncorrelated noises algorithm as p ma x increases. 



NUMBER OF SAMPLES L 

Fig. 6. Combining gain versus L with p = 0.7 for worst-case Rg- 


VI. Conclusion 

The correlated noises signal-combining scheme developed in this article has been shown to be asymp- 
totically optimal as the number of residual carrier-signal samples used in the estimates of the optimum 
weight coefficients increases. The numerical examples considered here show that convergence of the com- 
bining algorithm’s SNR performance to the optimum achievable SNR performance level to within 0.1 dB 
occurs at about L = 1000 signal samples. Hence, real-time operation is possible, although the inversion 
of a complex-valued A-dimensional matrix is required at every update of the weight coefficient estimates. 
The numerical examples considered here consistently demonstrate a significant performance superiority 
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GAIN IN SNR, dB GAIN IN SNR, dB 



FRACTION OF TOTAL POWER 


Fig. 7. Combining gain versus P with L = 5000 for worst-case Bb- 



—2.5 I t i t i t i i i i I i i i « I i i i i I 

0.00 0.05 0.10 0.15 0.20 


MAXIMUM CORRELATION COEFFICIENT 
Fig. 8. Combining gain versus p max with L = 5000 for worst-case Rg. 
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of the correlated noises combining scheme over the uncorrelated noises signal-combining algorithm. It 
appears to be the combining system of choice in the presence of significant correlations between the noise 
processes in different array feed received signals. A degree of caution should be exercised in extrapolating 
the expected amount of SNR gain or loss from the numerical examples considered in Section V. The best- 
case and worst-case performances in these examples should not be viewed as being typical. Moreover, 
these two cases also should not be viewed as being best- and worst-case performances in general. 
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Appendix A 

Combiner Output Signa! Variance 


An explicit expression for Var [s c (i B )] is derived in this appendix. For convenience in the following 
analysis, let us use W M I as the shortened notation for W m l Ua), X_ml f° r 2Lml(}a), and R A ML for 
Ba,ml(^a)- Using Eqs. (26) and (22), we can write 


Var[s c (i B )] =Var [w T ML VeP a ^ s ' 


=e [(wl L r) (k t ml v)] - | W%V 


=Z f E [W ML W T ML 


V - 1 wj>v | 2 


ssV^E 


(L-K- l) 2 --i * *-1 

cos 2 # -—aml2Lml*mlBa,ml 


v-\wlv\ 2 


(A-l) 


Since R a ,ml independent of X_ ML , we can write 


E 


• - -i - - 1 ~ -i 

Ba,ml2Lml2LmlBa,ml 


' =E [E [r a i 


ml^ml£-mlBa,ml I Ba,ml 




-E R Ay ML E I B A ,ML Ba,ML 


— E Sa.mlE [A ML A ML j R Ay ML 


(A-2) 


Moreover, since X ml Ua) is a complex Gaussian random vector with mean X_ and covariance matrix 
(1 /L)R a , we have 


E [x ML xi L ]=^R A + XX^ 


(A-3) 


Therefore, it follows from Eqs. (A-l), (A-2), and (A-3) that 


Var[s c (i B )] = 


{L-K- 1 )' 
L 2 77 2 cos 2 <5 


U f E 


Ra.ml (jBa+XX^R 


-i 

■A, ML 


V-\W%V\ 2 


(L-K-iy 

r] 2 cos 2 6 


U f E 


(jBa+ 2UL^ (LR AML y^V- | WlV | 2 (A-4) 
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Let Tr(A) denote the trace of a square matrix A and let I_k denote the K x K identity matrix. Recall 
that LR a ml is a K x K CW(R A , L — 1) distributed random matrix. It has been shown in [8] that if A 
is a K x K CW(R a , L — 1) random matrix and C_ is any constant K x K matrix, then 


E [- 1 (L — K)(L - K -2) 


Ea'CRA 1 + 


(L - K)(L -K- 1 )(L - K - 2) 


Tr^ 1 )^ 1 


(A-5) 


for L > K. So, using Eq. (A-5) in Eq. (A-4) gets 

V ” MW1 -I L-W-K-] " *> - *1 + *?**&) 

+ Tr (±l K + XX' Ea 1 ) Sa } Y. - | W£v | 2 

Ui [i -) +*&**&£ 

+ (L-K)iAl~-L^A ,Tr (r& + “’«') ~ ' I s (a-6) 


Since A = V. cos 6, and R s = r]R A , it follows from Eq. (29) that 


^ 2 - 6 (y}R A l X2C}R A l v) = I V'R-^V I 2 = I WZ,V I 2 = 7 2 MAX (A-7) 


Since R. A l is Hermitian, Tr( AA t R 1 1 ) = XJ R A X_- So we have, by using Eq. (A-7), 




?7 2 cos 2 (5 


-r+X^RT/X 


V'Ea'V 


K 

Lrf 2 cos 2 8 


(V-'Ea'y) +v ] Eb 1 Yv}Eb 1 v 


X 2 

f7 max + 7max 


Lrjcos 2 8 

Finally, using Eqs. (A-7) and (A-8) in Eq. (A-6) yields the following expression for Var[s c (i B )]: 


(A-8) 


Var [s c (i B )] 


(L - K - 1 )(L - 1) / \ 


+ 


(L - A)(L - K - 2) \Lr]cos 2 8 ) L — K — 2 


„,2 

A MAX 


(A-9) 


Note that Var[s c (i B )] approaches zero in the limit as L tends to infinity. This is because the estimate 
W_ ml {ia) converges with probability one to the optimum weight vector W Q given by Eq. (23) as L tends 
to infinity. 


50 



Appendix B 

Combiner Output Noise Variance 


Consider next an explicit expression for Var[n c (is)]. We shall also employ here the shortened notations 
W ML' X Ml • ns in Appendix A. Moreover, we shall use the shortened notation n for n(*s). Since 

W ML is independent of n, a derivation similar to that establishing Eq. (A-2) can be used along with 
Eq. (22) to get 


Var[n c (* B )] =E [wl L nrSW ML 


=E [W T ML R B W* ML ] 


-J2 ~ TJCOS 2 6^ aRa,MlX-ML (B-l) 


Since R a ,ml ls independent of X.ml> using the same approach on the expected value in Eq. (B-l) gets 




K' ml E 


{^B^aml) —a {^Ra,ml) 


%-ML 


(B-2) 


Next, using the property of Eq. (A-5) for the KxK CW(R A , L — 1) distributed random matrix LR aml 
in Eq. (B-2) results in the following expression: 


Var [n c (i B )\ 


(L-K- l ) ra r*t [(L-1)-K\E?+KEa* 

~ 77cos 2 <5 b [~ ML ( L - K)(L -K-2) ~ ML , 


(L-K-1){L-1) 
>cos 2 <5(L - K){L - K - 2) 


E 


[^mlBTaILml 


It follows from Eq. (A-3) that 

E |^A ML; R a 1 A ml | =E jTr [k a X.ml^Lml}^ 
=Tr (r^E [KmlKml]) 

=j+x 1 Ra 1 2L 


(B-3) 


(B-4) 
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Finally, using Eq. (B-4) in Eqs. (B-3) and (13) yields 


Var [n c (i B )\ 


1 + 


LK - K 2 - K + 1 


(L - K)(L - K - 2) J 


7 MAX + 


(L-1)(L-K-1) 


K 


(L - K)(L - K - 2)t ? cos 2 «5 L 


(B-5) 


Appendix C 

Mean-Square Convergence of Estimated Weights 


We employ the shortened notations W_ ML , ILmli an< ^ Qa,ml> similar to the usage in the previous 
appendices. Since 


W 


ML 


W r 


= Tr 


(Kml 


Ko) (] W ML 


w, 


»)1 


and 


E 




= E 


K ML W 


* 1 

ML 


- Wn Wl 


we need only establish that E \W ML w\ iL ] — > W 0 W p as L — > oo to prove the mean-square convergence 
of W.ml to Wo- Using Eqs. (17) and (18), we can write 


E 



(L - 1 - K) 2 
L 2 r] 2 cos 2 6 



-i . * -t 
X-mlZml 



It then follows from using Eqs. (A-2) and (A-3) in Eq. (C-l) that 


E 


W 


.mlW-ml 


{L-l- K)‘ 
; 2 6 


E 


Tj COS 


(lR^^ (^R* a + X*X T ^ [lR\ m ^ 1 


(C-l) 


(C-2) 


So, by using the property of Eq. (A-5) for the CW(R A ,L - 1) distributed random matrix LR aml in 
Eq. (C-2), we can write 


E 


\w ML wl L 


L-l-K 

T) 2 COS 2 8 


(L-l-K) 


\ mr 1 + mr l 2rx T (e^)- 1 ] + [f + * T (Ea)- 1 ^] mr 


(L - K)(L - K - 2) 


Therefore, E [W ML w} ML \ — > (RY)~ l X* ^(R^Y 1 lY cos 2 8 = W 0 W} 0 as L — > oo, thereby establish- 
ing the mean-square convergence of W MI to W 0 . 
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